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The semiclassical Einstein-Langevin equations which describe the dynamics of stochastic pertur- 
bations of the metric induced by quantum stress-energy fluctuations of matter flelds in a given state 
are considered on the background of the ground state of semiclassical gravity, namely, Minkowski 
spacetime and a scalar field in its vacuum state. The relevant equations are explicitly derived for 
massless and massive fields arbitrarily coupled to the curvature. In doing so, some semiclassical 
results, such as the expectation value of the stress-energy tensor to linear order in the metric per- 
turbations and particle creation effects, are obtained. We then solve the equations and compute the 
two-point correlation functions for the linearized Einstein tensor and for the metric perturbations. 
In the conformal field case, explicit results are obtained. These results hint that gravitational fluctu- 
, ations in stochastic semiclassical gravity have a "non-perturbative" behavior in some characteristic 

■ correlation lengths. 
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QQ ! I. INTRODUCTION 

<N : 

It has been pointed out that the semiclassical theory of gravity cannot provide a correct description of the 
\ dynamics of the gravitational field in situations where the quantum stress-energy fluctuations are important |j . 

QQ In such situations, these fluctuations may have relevant back-reaction effects in the form of induced gravitational 
fluctuations which, in a certain regime, are expected to be described as classical stochastic fluctuations. A 
. generalization of the semiclassical theory is thus necessary to account for these effects. In two previous papers, 
T-H ' Refs. and we have shown how a stochastic semiclassical theory of gravity can be formulated to improve the 
, description of the gravitational flcld when stress-energy fluctuations are relevant. 

In Ref . 1^ , we adopted an axiomatic approach to construct a perturbative generalization of semiclassical gravity 
which incorporates the back reaction of the lowest order stress-energy fluctuations in the form of a stochastic correction. 
O We started noting that, for a given solution of semiclassical gravity, the lowest order matter stress-energy fluctuations 
^ '■ can be associated to a classical stochastic tensor. We then sought a consistent equation in which this stochastic tensor 
5-H I was the source of linear perturbations of the semiclassical metric. The equation obtained is the so-called semiclassical 
Einstein-Langevin equation. 

J> , In Ref. 10 , we followed the idea, first proposed by Hu in the context of back reaction in semiclassical gravity, of 
' viewing the metric field as the "system" of interest and the matter fields (modeled in that paper by a single scalar field) 
rN \ as being part of its "environment." We then showed that the semiclassical Einstein-Langevin equation introduced in 

■ Ref. flCt] can be formally derived by a method based on the influence functional of Feynman and Vernon (see also 
■ ■ ■ Ref. |13|). That derivation shed light into the physical meaning of the semiclassical Langevin-type equations around 



speciflc backgrounds previously obtained with the same functional approach [gd 23|, since the stochastic source term 
was shown to be closely linked to the matter stress-energy fluctuations. We also developed a method to compute the 
semiclassical Einstein-Langevin equation using dimensional regularization, which provides an alternative and more 
direct way of computing this equation with respect to previous calculations. 

This paper is intended to be a first application of the full stochastic semiclassical theory of gravity, where we 
evaluate the stochastic gravitational fluctuations in a Minkowski background. In order to do so, we flrst use the 
method developed in Ref. [l0|] to derive the semiclassical Einstein-Langevin equation around a class of trivial solutions 
of semiclassical gravity consisting of Minkowski spacetime and a linear real scalar field in its vacuum state, which 
may be considered the ground state of semiclassical gravity. Although the Minkowski vacuum is an eigenstate of the 
total four-momentum operator of a field in Minkowski spacetime, it is not an eigenstate of the stress-energy operator. 
Hence, even for these solutions of semiclassical gravity, for which the expectation value of the stress-energy operator 
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can always be chosen to be zero, the fluctuations of this operator are non- vanishing. This fact leads to consider the 
stochastic corrections to these solutions described by the semiclassical Einstein-Langevin equation. 

We then solve the Einstein-Langevin equation for the linearized Einstein tensor and compute the associated two- 
point correlation functions. Even though, in this case, we expect to have negligibly small values for these correlation 
functions at the domain of validity of the theory, i.e., for points separated by lengths larger than the Planck length, 
there are several reasons why we think that it is worth carrying out this calculation. 

On the one hand, these are, to our knowledge, the first solutions obtained to the full semiclassical Einstein- 
Langevin equation. We are only aware of analogous solutions to a "reduced" version of this equation inspired in a 
"mini-superspace" model |^ . There is also a previous attempt to obtain a solution to the Einstein-Langevin equation 
in Ref. jlj], but, there, the non-local terms in the Einstein-Langevin equation were neglected. 

The Einstein-Langevin equations computed in this paper are simple enough to be explicitly solved and, at least for 
the case of a conformal field, the expressions obtained for the correlation functions can be explicitly evaluated in terms 
of elementary functions. Thus, our calculation can serve as a testing ground for the solutions of the Einstein-Langevin 
equation in more complex situations of physical interest (for instance, for a Robertson- Walker background and a field 
in a thermal state). 

On the other hand, the results of this calculation, which confirm our expectations that gravitational fluctuations are 
negligible at length scales larger than the Planck length, can be considered as a flrst check that stochastic semiclassical 
gravity predicts reasonable results. 

In addition, we can extract conclusions on the possible qualitative behavior of the solutions to the Einstein-Langevin 
equation. Thus, it is interesting to note that the correlation functions are characterized by correlation lengths of the 
order of the Planck length; furthermore, such correlation lengths enter in a non-analytic way in the correlation 
functions. This kind of non-analytic behavior is actually quite common in the solutions to Langevin-type equations 
with dissipative terms and hints at the possibility that correlation functions for other solutions to the Einstein-Langevin 
equation are also non-analytic in their characteristic correlation lengths. 

The plan of the paper is the following. In Sec. ||, we give a brief overview of the method developed in Ref. |Tc| ] 
to compute the semiclassical Einstein-Langevin equation. We then consider the background solution s of semiclassical 



gravity consisting of a Minkowski spacetime and a real scalar field in the Minkowski vacuum. In Sec. Ill, we compute 



the kernels which appear in the Einstein-Langevin equation. In Sec. IV, we derive the Einstein-Langevin equation 
for metric perturbations around Minkowski spacetime. As a side result, we obtain some semiclassical results, which 
include the expectation value of the stress-energy tensor of a scalar field with arbitrary mass and arbitrary coupling 
parameter to linear order in the metric perturbations, and also some results concerning the production of particles 
by metric perturbations: the probability of particle creation and the number and energy of created particles. In 
Sec. we solve this equation for the components of the linearized Einstein tensor and compute the corresponding 
two-point correlation functions. For the case of a conformal field and spacelike separated points, explicit calculations 
show that the correlation functions are characterized by correlation lengths of the order of the Planck length. We 
conclude in Sec. VI with a discussion of our results. We also include some appendices with technical details used in 
the calculations. 

Throughout this paper we use the (-1- -I- +) sign conventions and the abstract index notation of Ref. |Q, and we 
work with units in which c ^ h ~ 1 . 



II. OVERVIEW 

In this section, we give a very brief summary of the main results of Refs. and which are relevant for 
the computations in the present paper. One starts with a solution of semiclassical gravity consisting of a globally 
hyperbolic spacetime (A^,5ab), a linear real scalar field quantized on it and some physically reasonable state for this 
field (we work in the Heisenberg picture) . According to the stochastic semiclassical theory of gravity [^,|lO| , quantum 
fiuctuations in the stress-energy tensor of matter induce stochastic linear perturbations hab to the semiclassical metric 
gab- The dynamics of these perturbations is described by a stochastic equation called the semiclassical Einstein- 
Langevin equation. 

Assuming that our semiclassical gravity solution allows the use of dimensional analytic continuation to define 
regularized matrix elements of the stress-energy "operator," we shall write the equations in dimensional regularization, 
that is, assuming an arbitrary dimension n of the spacetime. Using this regularization method, we use a notation in 
which a subindex n is attached to those quantities that have different physical dimensions from the corresponding 
physical quantities. The n-dimensional spacetime (Al, ffab) has to be a solution of the semiclassical Einstein equation 
in dimensional regularization 
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(2.1) 



where Gs, A^, as and (3b are bare coupling constants and Gab is the Einstein tensor. The tensors I?"^ and B'^'^ 
are obtained by functional derivation with respect to the metric of the action terms corresponding to the Lagrangian 
densities RabcdR'^'"^'^ — RabR"'^ and i?^, respectively, where Rabcd is the Riemann tensor, Rab is the Ricci tensor and 
R is the scalar curvature (see Ref. |l^ for the explicit expressions for the tensors D"'' and B""^). In the last equation, 
Tjj'' is the stress-energy "operator" in dimensional regularization and the expectation value is taken in some state 
for the scalar field in the n-dimensional spacetime. Writing the bare coupling constants in Eq. (2.1) as renormalized 
coupling constants plus some counterterms which absorb the ultraviolet divergencies of the right hand side, one can 
take the limit n— >4, which leads to the p hysi cal semiclassical Einstein equation. 

Assuming that gab is a solution of Eq. (2.1), the semiclassical Einstein-Langevin equation can be similarly written 
in dimensional regularization as 
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^[g + h] = ^,-'^-^^(^f^'^)[g + h] + 2^,-^-^^C, (2.2) 



where hab is a linear stochastic perturbation to gab-. 
stochastic tensor characterized by the correlators 

(^f(2:))c=0. 



and h''^ = g'"'g'"^h. 
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In this last equation, is a Gaussian 



(CWC(2/)>c=iVr^''(x,y) 



(2.3) 



where 8N-'>-''ix,y) = ({C(^), C(y)}> [.9] 



with if^ 



^ n r 
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here, 



( )c means statistical average and { , } denotes 
an anticommutator. As we pointed out in Ref. |10|, the noise kernel jv°feg^(x, y) is free of ultraviolet divergencies in 
the limit n— i-4. Therefore, in the semiclassical Einstein-Langevin equation (2.2^ one can perform exactly the same 



renormalization procedure as the one for the semiclassical Einstein equation (2.1), and Eq. (2.2) yields the physical 
semiclassical Einstein-Langevin equation in four spacetime dimensions. 

In Ref. 0, we used a method based on the CTP functional technique applied to a system-environment interaction, 
more specifically, on the influence action formalism of Feynman and Vernon, to obtain an explicit expression for the 
expansion of {T^^)[g + h] up to first order in hcd- In this way, we can write the Einstein-Langevin equation (2.2) in a 
more explicit form. This expansion involves the kernel H^j'"^'^{x,y) = H^'"^'^{x,y) + H^'"^'^ {x , y) , with 
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(rf (x)r„=^(y)))[g], ffrr(x,y) ^ {[f^^x), f,fiy)])[g] 



(2.4) 



where [ , ] means a commutator, and we use the symbol T* to denote that we have to time order the field operators 
$„ first and then to apply the derivative operators which appear in each ter m o f the product T'^^{x)T'^'^{y), where 
T"'' is the classical stress-energy tensor; see Ref. [|l^ for more details. In Eq. (2^), all the ultraviolet divergencies in 
the limit n ^ 4, which shall be removed by renormalization of the coupling constants, are in some terms containing 
{^"^{x)) and in H^J^'^'^ {x , y) , whereas the kernels N^^'''^{x,y) and H^^'^{x,y) are free of ultraviolet divergencies. These 
two last kernels can be related to the real and imaginary parts of (t'^{x) i^f{y)) by 
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K'^'^ix^y) = -Re{i'^'{x)i'^\y)) ^ 



<^'(^,y) = Jim(ir(3:)tr(y)) 



(2.5) 



We now consider the case in which we start with a vacuum state |0) for the field quantized in spacetime {A4, gab) - 
In this case, it was shown in Ref. jlOj that all the expectation values entering the Einstein-Langevin equation ( |2.2| ) 
can be written in terms of the Wightman and Feynman functions, defined as 



G+(x,2/) = (0|<l„(a;)$„(2/)|0)[<7], 



zG,„(x,y) = (0|T(|.„(a:)<I.„(y))|0)[5]. 



(2.6) 



For instance, we can write ($^(a;)) = iGp^^{x,x) — G'^(x,x). The expressions for the kernels, which shall be used in 
our calculations, can be found in Appendix 



A. Perturbations around Minkowski spacetime 



An interesting case to be analyzed in the framework of the semiclassical stochastic theory of gravity is that of a 
Minkowski spacetime solution of semiclassical gravity. The flat metric r}ab in a manifold M = IR"* (topologically) and 
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the usual Minkowski vacuum, denoted as |0), give the class of simplest solutions to the semiclassical Einstein equation 
[note that each possible value of the parameters (m^,^) leads to a different solution], the so called trivial solutions 
of semiclassical gravity [^Sj. In fact, we can always choose a renormalization scheme in which the renormalized 
expectation value (0| |0) [77] = 0. Thus, Minkowski spacetime (H .rjab) and the vacuum state |0) are a solution 
to the semiclassical Einstein equation with renormalized cosmological constant A = 0. The fact that the vacuum 
expectation value of the renormalized stress-energy operator in Minkowski spacetime should vanish was originally 
proposed by Wald [|j and it may be understood as a renormalization convention There are other possible 

renormalization prescriptions (see, for instance, Ref. p^) in which such vacuum expectation value is proportional to 
and this would determine the value of the cosmological constant A in the semiclassical equation. Of course, all 
these renormalization schemes give physically equivalent results: the total effective cosmological constant, i.e., the 
constant of proportionality in the sum of all the terms proportional to the metric in the semiclassical Einstein and 
Einstein-Langevin equations, has to be zero. 

Although the vacuum |0) is an eigenstate of the total four-momentum operator in Minkowski spacetime, this state 
is not an eigenstate of 2^[77]- Hence, even in these trivial solutions of semiclassical gravity, there are quantum 
fluctuations in the stress-energy tensor of matter and, as a result, the noise kernel does not vanish. This fact leads 
to consider the stochastic corrections to thi s class of trivial solutions of semiclassical gravity. Since, in this case. 



the Wightman and Feynman functions (2^), their values in the two-point coincidence limit, and the products of 
derivatives of two of such functions appearing in expressions ( |All ) and (^) (Appendix are known in dimensional 



regularization, we can compute the semiclassical Einstein-Langevin equation using the method outlined above. 

In order to perform the calculations, it is convenient to work in a global inertial coordinate system {x^} and in 
the associated basis, in which the components of the flat metric are simply 77^^ — diag(— 1, !,...,!). In Minkowski 
spacetime, the components of the classical stress-energy tensor functional reduce to 

T'^'ir], $] = d^'^d"^ - ^ rj'"' dP^dp<i> - i Ty^'^m^^.^ + ^ (^M^^n _ Qf^d") 4>^ (2.7) 

where □ = d^d^, and the formal expression for the components of the corresponding "operator" in dimensional 
regularization is 

Triv] = I a"-*™} + v^'^l (2.8) 

where D^"^ are the differential operators "Di^" = (^ — l/4)i]^^'^0^ — £,d!^d^ and $„(a;) is the field operator in the 
Heisenberg picture in an n-dimensional Minkowski spacetime, which satisfies the Klein-Gordon equation (□ — m^) = 
0. 

Notice, from ( |2.S| ), that the stress-energy tensor depends on the coupling parameter ^ of the scalar field to the 
scalar curvature even in the limit of a flat spacetime. Therefore, that tensor differs in general from the canonical 
stress-energy tensor in flat spacetime, which corresponds to the value ^ = 0. Nevertheless, it is easy to see that 
the n-momentum density components ^'ri'(^) iv] temporary use this notation to indicate the dependence on the 

parameter ^) and T^^^^_q^ [77] differ in a space divergence and, hence, dropping surface terms, they both yield the same 
n-momentum operator: 

P'^l : M : = I d-^x : f O';^^^^ [77] : , (2.9) 

where the integration is on a hypersurface = constant (P'^ is actually independent of the value of x'^) and we use the 
notation for coord inat es a;^ = (x^jX), i.e., x are space coordinates on each of the hypersurfaces x° = constant. The 
symbol : : in Eq. ( |2.9| ) means normal ordering of the creation and annihilation operators on the Fock space built on 
the Minkowski vacuum |0) (in n spacetime dimensions), which is an eigenstate with zero eigenvalue of the operators 



(2.9) 



The Wightman and Feynman functions (p.6D in Minkowski spacetime are well known: 

G+{x,y) = {0\^,,{x)K{y)\OM=^A+{x~y), 

a„(x, y) = -I (0| T(|.„(x)l>„(y)) |0>[77] - A,„(a; - y), (2.10) 



with 



4 



r 



dJ^k 



■ikx 



(2.11) 
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where fc^ = r/^^k'^k'^ and kx = rj^^^k^^x^ . Note that the derivatives of these functions satisfy d^A^{x~y) = 9^A+(x— y) 
and 9^A+(a; — y) — — 9^A+(x — y), and similarly for the Feynman propagator Aj?„(x — y). 



To write down the semiclassical Einstein equation ( p.l[ ) for this case, we need to compute the vacuum expectation 
value of the stress-energy operator components (^). Since, from ( 2.10 ), we have that (0|l>^(x)|0) = iA^^(O) = 
«A+(0), which is a constant (independent of x), we have simply 



|T'n^nO)W = \ (0| {9^l>„,5''$„} |0)M = -z{d^^d^ 



ki'k'' 



cC'k 



(27r)" fc2 + m2 - i€ 
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i/2 



(2.12) 

where the integrals in dimensional regularization have been computed in the s tan dard way (see Appendix ^) and 
where T{z) is the Euler's gamma function. The semiclassical Einstein equation ( |2.l|) , which now reduces to 

Ab 



(2.13) 



simply sets the value of the bare couphng constant A^/Gs- Note, from (2.12), that in order to have (0| TjJ |0)[?7] =0, 
the renormalized (and regularized) stress-energy tensor "operator" for a scalar field in Minkowski spacetime has to 
be defined as 



2 (47r)2 \^47r^2 

which corresponds to a renormalization of the cosmological constant 



Ab 
Gb 



A 



G TT n{n-2) 



where 



1 



(n-4) \47r^2 



n-4 



2 



47r^2 



+ 0(n-4). 



(2.14) 



(2.15) 



(2.16) 



being 7 the Euler's constant. In the case of a massl ess scalar field, m^^O, one simply has As/Gs = A/G. Introducing 
this renormalized coupling constant into Eq. ( 2.13| ), we can take the limit n^A. We find again that, for (IR^,77ab, |0)) 
to satisfy the semiclassical Einstein equation, we must take A = 0. 

We are now in the position to write down the Einstein-Langevin equations for the components /i^j/ of the stochastic 
metric perturbation in dimensional regularization. In our case, using (0|<&2 (a^)|o) = zAp„(0) and the explicit expression 
for Eq. (|2.2|) found in Rcf. llG] , we obtain that this equation reduces to 



SvrG, 



G'^>'"' + Ab ( h>'''--rif"'h 



{x) - - agD^^^t^^ix) - 2/3bS<i>'"^(x) 



eG'^'^''(a;)Ai-("-^) zA,„(0) + 2 / d"y ^-("-4)iJ^-"/3(x, y) h^p{y) = 2^^''(x), 



(2.17) 

where ^'''^ are the components of a Gaussian stochastic tensor of zero average and 

{r{x)C''{y))c = fi-^^''-'^Nr"^{x,y), (2.18) 

and where indices are raised in /i^^ with the fiat metric and h = h^. We use a superindex (i) to denote the components 
of a tensor linearized around the flat metric . In the last expressions, N^'^°'^{x^ y) and H!^^"-^{x^ y) are the components 
of the kernels defined above. In Eq. ( 2.17 ), we have made use of the explicit expression for G'^'^"^. This expression 
and those for D*^'^'' and B'-^'"-^'^ are given in Appendix^; the last two can also be written as 



where T!^" is the differential operator T!^" = -q^^^x — d^d!^- 



(2.19) 
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III. THE KERNELS FOR A MINKOWSKI BACKGROUND 



The kernels Nl^'^"'^{x,y) and Hl^'^°'^{x,y) — H^^"'^ {x,y) + H^^"' {x,y) can now be computed using (p^) and the 



expressions (Al) and (A3). In Ref. we have shown that the kernel Hj^^'^"^ {x,y) plays the role of a dissipation 
kernel, since it is related to the noise kernel, Nj^'^"^ {x , y) , by a fluctuation-dissipation relation. From the definitions 



(2.4) and the fact that the Minkowski vacuum |0) is an eigenstate of the operator P^, given by (2.9), these kernels 
satisfy 



(3.1) 



A. The noise and dissipation kernels 



Since the two kernels (2^) are free of ultraviolet divergencies in the limit n^4, we can deal directly with 

M'^-^^^ix - y) EE hm (0| i'^^ix) if{y) |0)M. (3.2) 

n—*4: 

The kernels 4iV'^''°''(a;,y) = Re Mt""^^ {x-y) and 4i7i"'"'^(a;, y) =lmM^"''^^{x-y) are actually the components of the 
"physical" noise and dissipation kernels that will appear in the Einstein-Langevin equations once the re norm alization 
procedure has been carried out. Note that, in the renormalization scheme in which T^''[r]] is given by (2.14), we can 
write M^""^l^{x-y) = {Q\f'^''{x)f'^^{y) where the limit n — > 4 is understood. This kernel can be expressed in 

terms of the Wightman function in four spacetime dimensions. 



A+(a;) 



-2m 



(2^ 



e*'=^(5(A:2 + m2)0(fc°), 



(3.3) 



in the following way: 



M^'''"''^{x) = -2 [a^a("A+(x) ^'^^^''l\+{x)+V^"'{^°'l\+{x) d>^A+{x)) 
+ V°''^{^^'A+{x)^''A+{x)) +V^"'V°'>^{A+^{x))] . 



(3.4) 



The different terms in Eq. (3.4) can be easily computed using the integrals 



lip) 



(2^)' 



d{k^ + m^) ei-k°) S[{k - pf + m2] 0(fc° - /), 



1p)^ 



d^k 
(2^)4 



fc^i • • • fc'''- (5(r + m') 6{-k'') 5[{k - pY + m'] 6'(fc" - p'^), 



(3.5) 



with r — 1,2, 3, 4, given in Appendix^; all of them can be expressed in terms of I{p). We obtain expressions (C1)-(|C 
It is convenient to separate I{p) in its even and odd parts with respect to the variables p'^ as 



where Isi-p) = hip) and Ia{~p) = 

Up) 
Up) 



/(p)=/s(p)+/a(p), 
'Ia{p)- These two functions are explicitly given by 



(3.6) 



1 



8(27r)3 
-1 



e{-p'-Am')Jl + i^, 



signp 







V-4m^) Wl + 4y . 



8(27r)3 

Using the results of Appendix O, we obtain expressions (|C4[)-(|C6|) and, after some calculations, we find 



(3.7) 
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1 + 4— lip) 
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(2^ 



(3.8) 



where = S, — 1/6. The real and imaginary parts of the last expression, which yield the noise and dissipation kernels, 
are easily recognized as the terms containing Is{p) and Ij^{p), respectively. To write them explicitly, it is useful to 
introduce the new kernels 



2\ — 



1 



19207r J (27r 



iVB(a:;m2,A0 



1 



d^p 
2887r y (27r)4 



(V-4m2)Wl + 4^ 



m 

3A^+ — 



i?B(a;;m^AO 



19207r y (27r)4 



sign/0(V-4m2)Jl + 4i^ 1 + 4^ , 



e*^"^ signp" 



-p^ - 4m^) J 1 + 4 



3Ae 



(3.9) 



28877 7 (27r)4 
and wc finally get 

N^^'^^x, y) = i {iT^^^'-Tl^'' - .^r-^"^) ^A(a;-y; m") + .Fr^,"'5iVB(a;-y; A^, 

iir"''(^,y) - i(3.F,^(".F«^-.F^^:^)i?A(a:-y;m2)+.F^^:''i^B(a;-y;™^A0. (3.10) 

Notice that the noise and dissipation kernels defined in ( |3.9| ) are actually real because, for the noise kernels, only the 
cos pa; terms of the exponentials e*''^ contribute to the integrals, and, for the dissipation kernels, the only contribution 
of such exponentials comes from the isinpx terms. 

We can now evaluate the contribution of the dissip ation kernel components Hj^'^'^^ {x , y) t o the Einstein-Langevin 
equations (2.17) [after taking the limit n— >4]. From (3.1C), integrating by parts, and using (2.19) and the fact that, 
in four spacetime dimensions, D^'-''^'^{x) = (3/2) A'^"'''(x) [the tensor A"'' is obtained from the derivative with respect 
to the metric of an action term corresponding to the Lagrangian density CabcdC"'^'^'^, where Cabcd is the Weyl tensor, 
see Ref. for details], it is easy to see that 



2 / d^y Hr^^ix, y) h^p{y) = / d^y [Dj,{x~y- m^)A''>^'' {y) + D^{x-y- m\ A0S<^'^''(2/)] 



(3.11) 

These non-local terms in the semiclassical Einstein-Langevin equations can actually be identified as being part of 



B. The kernel H^°''^{x,y) 



The evaluation of the kernel components H^^"^ {x,y) is a much more cumbersome task. Since these quantities 
contain divergencies in the limit n — > 4, we shall compute them using dime nsional regularization. Using Eq. ( |A3| ), 
these components can be written in terms of the Feynman propagator (2.11) as 



where 



Kt''"^'^{x) EE -fI''^''''^H2^^'^'^°'A^Jx)^'^^^''ApJx) + 2Vt"' (d'^A^Jx) d'^ A^Jx)) 



(3.12) 



+ 22?"^(5^A,„(x) 9''A,„(a;)) + 2V^^''V''^ {Al^{x)) 



j^t^^Qio^/s.^^i^x) 5'^' + ?7"''5(^A^„(a;) d"'^ 



- A,„(0) (r^^^V^^ + r^V^"") + 4 V^^r" (A.„(a;)n - m^A^jO)) 



8-{x) 



(3.13) 
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Let us define the integrals 

Up) ^ M"^""'^ 



1 



(27r)" (fc2 + m2 - ie) [{k - pY + - ie] 



/jMl . . . /cMr 



(27r)" (fc2 + to2 - ie) [{k - p)"^ + 



— ze 



(3.14) 



with r = 1, 2, 3, 4, and 



= „-("-4) 



1 



,M......^^-(n-4) 



dJ^k 

(27r)" (fc2 + to2 - ie) ' 



(27r)" (fc2 + _ ie) ' 



(3.15) 



with r=l,2, where a limit e^0+ is understood in all these expressions. Then, the diffe rent t erms in Eq. ( 3.13 ) can 
be computed using Eqs. (Dl)-(p^). The results for the expansions of the integrals (3.14) and (3.15) around n = A are 
given in A ppen dix |^. In fact, jg^ — and the remaining integrals can be w ritten in terms of /on Jnjp ) giv en in 
Eqs. (Bl) and (B4). Using the results of Appendix^, we obtain Eqs. (D7) and (D8) and, from Eqs. (D4 )-(□(:), we 
get 



-(«-4) 



(3.16) 



We are now in the position to work out the explicit expression for K^^°'^{x)^ defined in ( 3.13| ). We use 
Eqs. ( pi| ), the results (^, (M, (^ and (Q, the identities (5"(x) = (27r)-" /d>e*''"^, Jd^e'P'' f{p) = 
—J dPp e^P^ f{p) p^pt^" and d^d^jdJ^p e*^^ f{p) — — J (P^p e'^^ f{p) p^p'^ , where f{p) is an arbitrary function of and 
P'^'^ is the projector orthogonal to p^ defined as p^P'^'^ = ri'^'^p^ —p'^p", and the expansions in (Bl) and (B4) for Jn{p) 
and Jo„. After a rather long but straightforward calculation, we get, expanding around n = 4. 



2 



(47r)2 



3 (n-2) 

4m'^ 
n(n--2) 



{2'q^'(°'r)^>-r)^"'r^°'l^)5''{x) 



1 on a; ■'a; x x ) 



180 



(27r)" 

^ (3.F,^(".Ff)'' - TH'^Tf) + ^ (60^-11) .Fr^r'^' 



3AC+^ </>(p2) 



2 1 



6-ix) 
A^ix) }+0{n~A), 



where k„ and (pip^) have been defined in ( 2.16| ) and (B9), and A„(a;) is given by 



A„(x) 



(27r)" 



(3.17) 



(3.18) 



The imaginary part of ( ^.17| ) [which, using ( 3.12 ), gives the kernel components H^J^"^ {x,y)] can be easily 

obtained multiplying this expression by —i and retaining only the real part, Lpip^), of the function ipiP^)- Making use 
of this result, it is easy to compute the contribution of these kernel components to the Einstein-Langevin equations 
Integrating by parts, using Eqs. (E1)-(E5) and Eq. (2.19), and taking into account that, from Eqs. (2.12) and 



( 2.17 ) 
(2.13) 



m 



BttG, 



47r2 n{n-2) 



K„ + 0{n - 4), 



(3.19) 
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we finally find 

1 



A 



B 



90 



(27r)" 



Jp{x-y) 



1 



2 ' 



16 



28807r2 15 



(47r)2 
'1 



2 



6 



3 (n-2) 
10AnS<^>^''(x) 



1 + 4— i?'i)^-(y) + lO 3A^+— B<^»"'(y) 



A„(x - y) (8i?<^'^''(y) + SB'^'^^ly)) I + 0(n - 4). 



(3.20) 



C. Fluctuation-dissipation relation 



From expressions (3.10) and ( |3.9| ) it is easy to check that there exists a relation between the noise and dissipation 
kernels in the form of a fluctuation-dissipation relation which was derived in Ref. |jl^ in a more general context. 
Introducing the Fourier transforms in the time coordinates of these kernels as 



N''''"l^{x,y) 



2n 



and similarly for the dissipation kernel, this relation can be written as 

^M-/3(pO. y) ^ signpO 7V^''"^(pO; x, y), 

or, equivalently, as 

i/r"^(xO,x;2/0,y) = -- T dz^ iV^''"'3(z°, x; yO, y), 

where P( l/a: ") denotes the principal value distribution. 

From (3.1), taking the limit n^4, we see that the noise and dissipation kernels must satisfy 



J d^xN°^'°'^{x,y) = J d^xH°^'°''^{x,y) = 0. 



(3.21) 



(3.22) 



(3.23) 



(3.24) 



In order to check the last relations, it is useful to write the J-!^'^ derivatives in expressions (3.10) using 
^x'^ I d'^pe^P^^~y^ f{p) — — J d'^pe^P^^~y^ f {p) p^ P'^" , where f{p) is any function of and P'^'^ is the projector orthog- 
onal to p^ defined above. The identities ( |3.24 ) follow by noting that p^P^^ = —p^Pi and p'^P'^^ — —p^p"^, where we 
use the index i = 1, 2, 3 to denote the space components, and that / d'^x exp(ipia;*) = (2Tr)^ Il^^i Hp^)- It is also easy 
to check that the noise kernel satisfies d^N'^'^°'^{x,y) — and, hence, the stochastic source in the Einstein-Langevin 
equations will be conserved up to first order in perturbation theory. 



IV. THE SEMICLASSICAL EINSTEIN-LANGEVIN EQUATIONS 



The results of the previous secti on are now ready to be introduced i nto the Einstein-Langevin equations (2.17). In 
fact, substituting expression ( 3.20 ) in such equations, and using Eqs. (D4) and (Bl) for the /i"*^""*^ Af„(0) term, we 
get 



^ G"'-(.) - I asD-^^^{x) - 2/3^B<^>-(.) + ^ 



m 1 

^(n-2) ^90 



{x) 



J --]^D'■'^^"'(x) +(- -10A^^S<^)^'^(a;) 

28807r2 I 15 ^ ' \& ^) ^ ' 
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(27r)" 



Jp{x-y) 



2 \ ^ / 2 \ ^ 



m 



d"yA„(a;-y) (87^<^>^'^ + 5B<^>^'')(2/) I + 2 / d^y fi^'^^-^^ Hj^^;"^ {x,y) /^^.^(y) + 0(n-4) 



(4.1) 



Notice that the terms containing the bare cosmological constant have canceled. These equations can now be renor- 
malized, that is, we can now write the bare couphng constants as renormahzed coupUng constants plus some suitably 
chosen counterterms and take the limit n^4. In order to carry out such a procedure, it is convenient to distinguish 
between massive and massless scalar fields. We shall evaluate these two cases in different subsections. 



A. Massive field (m / 0) 



In the case of a scalar field with mass m 0, we can use, as we have done in Eq. ( 2.15| ) for the cosmological 
constant, a renormalization scheme consisting on the subtraction of terms proportional to k„. More specifically, we 
may introduce the renormahzed coupling constants 1/G, a and /3 as 



1 _ 1 

Gb G 

as = a + 
I3b = P + 



TT ^ (n-2) 
1 1 



K„ + 0(n-4). 



— K„ + 0(ji-4), 
2 K„ + 0(n-4). 



(47r)2 120 
32^ 



(4.2) 



Note that for conformal coupling, = 0, one has 1/Gb = 1/G and Pb ~ that is, only the co uplin g constant a 
and the cosmological constant need renormalization. Substituting the above expressions into Eq. ( [4.l[ ), we can now 
take the limit n-^4, using Eqs. ( |3.1^ ), (|3.1l| ) and the fact that, for n = 4, /^(^""'(a;) = (3/2) ^(^'^''(x). We obtain 
the semiclassical Einstein-Langevin equations for the physical stochastic perturbations h^i, in the four-dimensional 
manifold A^=IR^. Introducing the two new kernels 



19207r2 J (27r)4 



-jTTSignp 







V-4m2)Wl + 4 — 



H^{x-m\AO 



1 



2887r2 J (27r)4 



3AC + -:t 



3 p2 y 

iTTsignp° 0{-p^-'im^) 



1 + 4^ 



1 

6 p2 



(4.3) 



where (pip^) is given by the restriction to 7i = 4 of expression ( B1C| ), these Einstein-Langevin equations can be written 



as 



SttG 



&'^'"'{x) - 2 (aA<^'^'^(x) + /?S<^'^''(a:)) + 



1 



28807r2 



~-A'^'^^"'{x)+(--lQAi\B^''^"'{x) 



+ jd''y[HAix-y;m^)A'^'^^''{y) + HB{x~y;m^A0B^'^^^''{y)]^2e''ix), (4.4) 

where ^'^'^ are the components of a Gaussian stoc hastic tensor of vanishing mean value and two -poi nt correlation 
function {^'"'{x)C°''^{y))c = Nf"'"f^{x,y), given in ( |3.10| ). Note that the two kernels defined in ( [4.3| ) are real and 
can be split into an even part and an odd part with respect to the variables x^^, with the odd terms being the 
dissipation kernels Dx{x;m?) and DB{x;m?,A£^) defined in ([3.9[). In spite of appearances, one can show that the 
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Fourier transforms of the even parts of these kernels are finite in the hmit p^— s-0 and, hence, the kernels Hx and 
are well defined distributions. 

We should mention that, in a previous work in Ref. jist , the same Einstein-Langevin equations were calculated 
using rath er d ifferent methods. The way in which the result is written makes difficult a direct comparison with our 
equations (4.4). For instance, it is not obvious that in those previously derived equations there is som e a nalog of the 
dissipation kernels related to the noise kernels by a fluctuation-dissipation relation of the form (3.22) or (3.23). 



B. Massless field (m = 0) 



In this subsection, we consider the limit to— s-O of equations (4.1). The renormalization scheme used in the previous 
subsection becomes singular in the massless limit because the expressions (4.2) for as and /3_b diverge whe n to— > 0. 
Therefore, a different renormalization scheme is needed in this case. First, note that we may separate k„ in ( 2.16| ) as 
= + ^ ln(TO^//i^) + 0(n — 4), where 



(n-4) \47r 



7 \ K— 4 



.-4 2 \AtzJ 



+ 0(n-4). 



and that [see Eq. ( plO| )] 



lim [ip[p^) + ln(TO2//?)] = -2 + In 



(4.5) 



(4.6) 



Hence, in the massless limit, equations (4.1) reduce to 



1 



1 



SttG; 



■ &'^^^{x) - - aBD^'^^^ix) - 2/3bS<^»"'(x) + (^„-l) 



(47r)2 



90 



J ^i^zjCDf'^r-j;) + (^1 -10Ae^S<^'^''(x) + [d"y f 

28807r2 I 15 ^ ' \6 J J 

■90A^2^<^>'^''(y) 



(27r)" 



e'P^^'^y^ In 
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(x) 



These equations can be renormalized by introducing the renormalized coupling constants 1/G, a and /3 as 



1 1 



^ ^(^„-l) + 0(n-4), /35=;3+^(^„-l) + 0(71-4). 



(47r)2 120 



(4.7) 



(4.E 



Thus, in the massless limit, the Newtonian gravitational constant is not renormaliz ed a nd, in the conformal coupling 
case, = 0, we have again that Pb=P- Introducing the last expressions into Eq. (4.7), we can take the limit n^A. 
Note that, by making to = in (|3.9D, the noise and dissipation kernels can be written as 



NA{x\m^^Q) ^ N{x), 
DA{x;m^^O) = D{x), 



NB{x;m^^O,A^) = 6QA^^N{x), 
DB{x;m^:^0,A^) = 60A^2^(a;), 



where 



N{x) 



1 



19207r J (27r)4 
It is now convenient to introduce the new kernel 



D{x) 



d-^p 

19207r J (27r)4 



e'P^ signp" 



(4.9) 



(4.10) 



H{x■,^l') 



2N — 



1 



d-^p 

19207r2 J (27r)4 



lim 



19207r2 .-.0+ J (27r) 



d*p 



In 



e*P^ In 



— ITT sign p° 6{—p^) 
(p° + ie)2 



(4.11) 



Again, this kernel is real and can be written as the sum of an even part and an odd part in the variables x^, where 
the odd part is the dissipation kernel D{x). The Fourier transforms (4.1C) and (4.11) can actually be computed and. 
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thus, in this case, we have exphcit expressions for the kernels in position space. For N{x) and D{x), we get (see, for 
instance, Ref. P% ) 



N{x) 



1 



19207r 



D{x) 



1 



19207r3 



sign a; 



,0 



■S{x^). 



(4.12) 



where Vf denotes a distribution generated by the Hadamard finite part of a divergent integral (se e Refs. |g8| for the 
definition of these distributions). The expression for the kernel H{x\ /Lt^) can be found in Refs. |2^,0 and it is given 
by 



i/(a;;/i2) 



960 
1 



)^^6{x^)]+{l-^-\nn)5\x) 



960.2 \ I ^(-°) ^(Ixh A) ^ 8ix^) + [1-, - ln(,A)] 8\x) 



(4.13) 



See Ref. [g9| for the details on how this last distribution acts on a test function. Finally, the semiclassical Einstein- 
Langevin equations for the physical stochastic perturbations h^^ in the massless case are 



1 



■ G^^^^'''{x) - 2 (aA<i>^'^(x) + PB^^^f'" {x)) + 



1 



SttG ' ' ^ ' • ^ " • 28807r2 

+ [d'yH{x-y-iJi^) [A^'^^^^{y)+m^eB'''^''{v)\ = 2C^''(x), 



^v4(^'^''(a;)+(i-10AnS<^>''''(a;) 



(4.14) 



where the Gaussian stochastic source components have zero mean value and 



{i^''{x)e^{y))c^ lim iV'^"^"'5(x,2;) 

m— *0 



N{x-y). 



(4.15) 



It is interesting to consider the conformally coupled scalar field, i.e., the case A^ = 0, of particular interest because 
of its similarities with the electromagnetic field. It was shown in Refs. P,p^ that, for this field, the stochastic source 
tensor must be "traceless" (up to first order in perturbation theory around semiclassical gravity), in the sense that 
the stochastic variable ^'^ = 77^,^^^'^ behaves deterministically as a vanishing scalar field. This can be easily checked by 
noticing, from Eq. (|l|), that, when A^ = 0, one has {^jiix)^"'^ {y))c = 0, since :^/;' = 3n and JP'^^JP/^ = □JP"'^. The 
semiclassical Einstein-Langevin equations for this particular case [and generalized to a spatially flat Robertson- Walker 
(RW) background] were first obtained in Ref. (in this reference, the coupling constant /3 was set to zero). In order 
to compare with this previous result, it is worth noticing that the description of the stochastic source in terms of a 
symmetric and "traceless" tensor, with nine independent components ^'^'^ , is equivalent to a description in terms of a 
Gaussian stochastic tensor with the same symmetry properties as the Weyl tensor, with components ■C^''"'', defined as 
= —2dadi3^j^°"'^] this tensor is used in Ref. The symmetry properties of the ^a"'"/? ensure that there are also 
nine inde pend ent components in — 29Q,9^f^'"'''^. It is easy to show that, for this combination to satisfy the correlation 
(|4.15| ) with A^ = 0, the relevant correlators for the new stochastic tensor must be 



rclat 



(er"'^(a;)Cr^''(2/)>ec = T^^'^^P'^^^Nix - y), 



(4.16) 



where j't"^'^Pp<^>-^ jg g, linear combination of terms like r]'^'-Pri'^'^r]°'^r]/^^ in such a way that it has the same symmetries 
as the product of two Weyl tensor components (jt^^<^0 (j'"^^^ , its explicit expression is given in Ref. |jl^. Thus, after 
a redefinition of the arbitrary mass scale ^ in Eq. ( 4.14 ) to absorb the constants of proportionality of the local terms 
with A'^'^'^(x), one can see that the resulting equations for the A^ = case are actually equivalent to those found in 
Ref. M. 



C. Expectation value of the stress-energy tensor 



From the above equations one may extract the expectation value of the renormalized stress-energy tensor for a scalar 
field in a spacetime (IR^, r/ab -t- hab), computed up to first order in perturbation theory around the tri vial s olution of 
semiclassical gravity. Such an expectation value can be obtained by identification of Eqs. (4^) and ( 4.14 ) with the 
components of the physical Einstein-Langevin equation, which in our particular case simply reads 



12 



1 



SttG 



G 



(1)^11/ 



2(aA 



T^'')[v + h] + 2e''- 



(4.17) 



By comparison of Eqs. (4.4) and ( 4.14 ) with the last equation, we can identify 

1 



T^''ix))[r^ + h] = 



28807r2 



■ A<^"^''(a;)-Q-10Ae^S'^''"'(a;) 



(4.18) 



for a massive scalar field, m^O, and 



T^''ix))[v + h] = 



1 



5 V 6 



28807r2 

d^yH{x-y-i?) [A<^''"'(y) + 60A^25<^''"'(y)] + 0(/i2 



(4.19) 



for a massless scalar field, m = 0. Notice that in the massive case we have chosen, as usual, a renormalization scheme 
such that the expectation value of the renormalized stress-energy tensor does not have local terms proportional to the 
metric and the Einstein tensor The result (4.19) agrees with the general form found by Horowitz | ]30| , |3l[ using an 
axiomatic approach and coincides with that given in Ref. psf . The particular cases of conformal coupling, A^ = 0, 
and minimal coupling, A^ = — 1/6, are also in agreement with the results for this cases given in Refs. ||30|-|34| (modulo 
local terms proportional to and _B<i>^'' due to different ch oices of the renormalization scheme). For the case of 

a massive minimally coupled scalar field, A^ = — 1/6, our result ( 4.18 ) is equivalent to that of Ref. ||35|| . 

As it was pointed out above, in the case of conformal coupling, both for massive and massless scalar fields, one 
has =/3- This means that, in these cases, the terms proportional to B'-^^'^" in the above expectation values of the 
stress-energy tensor are actually independent of the renormalization scheme chosen. Due to the conformal invariance 
of / d'^x yf—gCcabdC""''"^, the tensor A"''' is traceless and we have A'-^'"^ = 0. Therefore, the terms with B'^^^^^" are 
precisely those which give trace to the expectation value of the stress-energy tensor in (4.18) and ( 4.19| ). In the 
massless conformally coupled case, m = and A^ = 0, such terms give the trace anomaly 0| up to first order in hi^^: 



h] 



^ 1 ^(DP 

28807r2 6 ^ 



28807r2 



(4.20) 



where we have used expression 



for B<^"^''. 



D. Particle creation 



We can also use the result ( 3.1C ) for the noise kernel to evaluate the total probability of particle creation and the 
number of created particles for a real scalar field in a spacetime (]R^,7ya6 + ^afc)- The metric perturbation hah (here 
an arbitrary perturbation) is assumed to vanish, either in an exact way or "asymptotically," in the "remote past" and 
in the "far future," so that the scalar field has well defined "in" and "out" many particle states. In that case, the 
absolute value of the logarithm of the vacuum persistence probability |(0,out|0,in)|2, where |0,in) and |0,out) are, 
respectively, the "in" and "out" vacua in the Heisenberg picture, gives a measure of the total probability of particle 
creation. On the other hand, the number of created particles can be defined as the expectation value in the "in" 
vacuum of the number operator for "out" particles. As it was shown in Ref. [^, the total probability of particle 
creation and one half of the number of created particles coincide to lowest non-trivial order in the metric perturbation, 
these are 



P[h] J c 



d^xd^y h^,{x) N'^-'^'^ix, y) Kp{y) + 0(;i^), 



(4.21) 



where N'^'^"^ {x,y) is the noise kernel given in ( |3.10D , which in the massless case reduces to ( [4.15| ). The above 
expression for the total probability of pair creation by metric perturbations about Minkowski spacetime was first 
derived in Ref. Using (KM, we can write P[h]=PA[h] + Psih] + 0{h^), where 



PA[h] = ^J d^xd^y {3TrK^ - Ti''Tf)NA{x-y;m'') h^,{x) Kpiy)^ 
P^[h] = J d^xdST^-TfNB{x-y-m\Ai) h^,{x) h^p{y). 



(4.22) 
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Integrating by parts (we always neglect surface terms), using expression (E5) for R^^\ which can also be written as 
^(1) = .jTM^'/j we find 



PB[h] - Jd^xd''yR'''>ix)NBix-y;m^,AOR''\y). 



(4.23) 



In order to work out PA[h], it is useful to take into account that, using the symmetry properties of the Weyl and 
Riemann tensors and the expression (E6) for one can write 



(i)paXT 



Cl'^,^{x)R^^^'"^'^{y) = -2C'^^'"^'^{x)S"^S'',d,drh^0{y). 



(4.24) 



Using the last identity, the expression (E7) for C'<^""^'^'^ and integrating by parts the first expression in ( 4.22 ) we get 

PaN= fd^xd^yCl'l^^{x)NA{x-y;m^)a'^^'^'^^{y). (4.25) 



Thus, PA[h] and Pslh] depend, respectively, on the Weyl tensor and the scalar curvature to first order in the metric 
pertu rbat ion. The result for the massless case, m — 0, can be easily obtained from the above expressions, using 
Eqs. (p^). If, in addition, we make = 0, i.e., conformal coupling, we have Pb[^] =0. Hence, for a conformal scalar 
field, particle creation is due to the breaking of conformal flatness in the spacetime, which implies a non-zero Weyl 
tensor. 

In order to compare with previously obtained results, it is useful to introduce the Fourier transform of a field ,f{x) 



as f{p) = Jd^a: e~'P^/(a;). Note that, if f{x) is real, then f{—p) = f*{p)- Using the expressions (3^) for the kernels 
Na and iVe, the above result for the total probability of particle creation and the number of particles created can 
also be written as 



P[h] 



1 



19207r J (27r 



d*p 



20 

y 



R''\p) 



3Ae + — 



Oih% 



(4.26) 



in agreement with the results of Ref. ||3^ (except for a sign in the coefficient of the term with |i?'^'(p)p). It is also 
easy to see that the above result is equivalent to that found in Ref. |Q if we take into account that, for integrals of 
the form / = / d'^p fai-.-a^ip) G{p'^) f*°'^'°"'{p), where fai---aA^) is any real tensor field in Minkowski spacetime and 
G{p^) is any scalar function of p^, one has that 



2 d^p9{p°)fa,...aAp)G{p')f' 



{p)^2 d^p0i^p")fa,...aAp)G{p')f 



ip)- 



(4.27) 



In the massless conformally coupled case, m — O and A^ = 0, the result (4.26) reduces to that found in Ref. [ p9[ . 

The energy of the created particles, E[h], defined as the expectation value of the "out" energy operator in the "in" 
vacu um c an be computed using the expressions derived in Ref. M- We find that this energy is given by an expression 
like ( 4.26 ), but with a factor 2p^ dip *^) in serted in the integrand |37|Jl0t] . Since the kernels Na and Da are rela ted b y 
the fluctuation-dissipation relation ( 3.22] ), and the same holds for Nb and Db, it is easy to see [similarly to ( 4.27 )] 
that 



Cl'L(3iP)G'''*""'^ip)DAip) 



R''\p) Db{p) 



0{h% 



(4.28) 



where Da{p) and Db{p) are the Fourier transforms of the dissipation kernels defined in ( |3.9[ ). For perturbations of a 
spatially flat RW spacetime (i.e., h^i, — 2lS.a{rf)ri^^, where x^ = r] is the conformal time and Aa(77) is the perturbation 
of the scale factor) , this last expression agrees with that of Ref. Q , see also Ref. . 

So far in this subsection the metric perturbations are arbitrary. We may also be interested in the particles created 
by the back reaction on the metri c du e to t he st ress-energy fluctuations. Then we would have to use the solutions of 
the Einstein-Langevin equations (4.4) and (4.14) in the above results. However, to be consistent, one should look for 
solutions whose moments vanish asymptotically in the "remote past" and in the "far future." These conditions are 
generally too strong, since they would break the time translation invariance in the correlation functions. In fact, the 
solutions that we find in the next section do not satisfy these conditions. 
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V. CORRELATION FUNCTIONS FOR GRAVITATIONAL PERTURBATIONS 



In this section, we solve the semiclassical Ein stein -Langevin equations (4^) and ( 4.14 ) for the components G*^''"' 
of the hnearized Einstein tensor. In subsection VA we use these solutions to compute the corresponding two-point 
correlation functions, which give a measure of the gravitational fluctuations predicted by the stochastic semiclassical 
theory of gravity in the present case. Since the linearized Einstein tensor is invariant under gauge transformations 
of the metric perturbations, these two-point correlation functions are also gauge invariant. Once we have computed 
the two-point cor relation functions for the linearized Einstein tensor, we find solutions for the metric perturbations 
in subsection V C and we show how the associated two-point correlation functions can be computed. This procedure 
to solve the Einstein-Langevin equations is similar to the one used by Horowitz js^ , see also Ref. , to analyze the 
stability of Minkows ki sp acctime in semiclassical gravity. 

From expressions (E2) and (E3) restricted to n — A, it is easy to see that A'^'^" and _B<i)^'' can be written in terms 
of G'^'^'' as 



3 



(5.1) 



where we have used that BD = T^. Therefore, the Einstein-Langevin equations (L4) and ( 4.14 ) can be seen as linear 
integro-differential stochastic equations for the components G'-^^^'^ . Such equations can be written in both cases, m^O 
and m = 0, as 



1 



G<^'^"(x) - 2 {aA'^^^"'{x) + pB''^^"'{x)) + / d^y [HA{x-y)A^'^^''{y) + HB{x-y)B^'^^''{y)] = 2C''{x), (5.2) 



SttG 

where the new constants a and and t he ke rnels Ha{x) and Hb{x) can be identified in each case by comparison 
of this last equation with Eqs. (L4) and ( 4.14| ). For instance, when m = 0, we have Ha{x) ~ H{x;iJ,'^) and Hb{x) = 
60A^^H{x; is^). In this case, we can use the arbitrariness of the mass scale n to eliminate one of the parameters a or 

In order to find solution s to these equations, it is convenient to Fourier transform them. Introducing Fourier 
transforms as in subsection IVD, one finds, from (5.1), 



i(i)A"'(p) = 2p2G<^'^''(p) - -p'P^"G<^'J^(p), 



^'^>'"'(p) = -2p2pM''G(i)2(p). 



Using these relations, the Fourier transform of Eq. (|^) reads 

F^';^(p)G<^'"'5(p) = 167rGe"'(p), 

where 

F^';^(p) ^ F,{p) S^Jf,) + F^{p)p'P^^^^0, 

with 



Fi(j3) = 1 + 167rGp2 FA(p)-2a 



F2{P) 
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ttG 



#A(p)+3FB(p)-2a-6^ 



In Eq. (5.4), i^^{p), the Fourier transform of ^^''(x), is a Gaussian stochastic source of zero average and 

where we have introduced the Fourier transform of the noise kernel. The explicit expression for N^^'^°'^{p) is found 
from ( |3l0| ) and (U) to be 



(5.3) 

(5.4) 
(5.5) 

(5.6) 
(5.7) 



N'"'°''^{p) 



1 



28807r 



(V-4m2) Jl 



10( 3Ae+ ) {p^)^pt^'^p'^l^ 



which in the massless case reduces to 

1 



lim N^"'°f^ 

m-tO 



ip) 



19207r 



di-p' 



(5.8) 



(5.9) 
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A. Correlation functions for the linearized Einstein tensor 



In general, we can wri te G*^'^"^ — (G'*^*^'')c + Gj-^'^'', where Gj-^"''^ is a solution to Eqs. ( pT^ ) [or, in the Fourier 
transformed version, (5^)] with zero average. T he a vera ges ( G'^'^"^)c must be a solution of the linearized semiclassical 



Einstein equations obtained by averaging Eqs. (5.2) 
case, m — 0) have been studied by several authors 



5.4)1. Solutions to these equations (specially in the massless 
2p^j34|,|25| , particularly in connection with the issue of 



the stability of the trivial solutions of semiclassical gravity. The two-point correlation functions for the linearized 
Einstein tensor are given by 



g^'"'^{x,x') = (G<^>^"(x)G<^'"''(x'))c - (G<^'^''(a;)),(G<^'"''(a;'))c = (G<'""'(x)G^''"^(x'))c 



(5.10) 



Next, we shall seek the family of solutions to the Einstein-Langevin equations which can be written as a linear 
functional of the stochastic source and whose Fourier transform, G'-^^^'^{p), depends locally on (."^{p). Each of such 
solutions is a Gaussian stochastic field and, thus, it can be completely characterized by the averages (G<^'^'')c and 
the two-point correlation functions ( ^.10 ). For such a family of solutions, Gf^"^''(p) is the most general solution to 
Eq. (5.4) which is linear, homogeneous and local in (."^{p). It can be written as 



(5.11) 



where D^'^^p[p) are the components of a Lorentz invariant tensor field distribution in Minkowski spacetime [by 
"Lorentz invariant" we mean invariant under the transformations of the orthochronous Lorentz subgroup; see Ref. |3C| ] 
for more details on the definition and properties of these tensor distributions] , symmetric under the interchanges a <-*■ /3 
and /x ^ J/, which is the most general solution of 



(5.12) 



In addition, we must impose the conservation condition to the solutions: p^Gf (p) — 0, where this zero must be 
understood as a stochastic variable which behaves detcrministically as a zero vector field. We can write D^'^ ^piP) ~ 
apip) + ^h" apiP)^ where Dp^^/iip) ^ particular solution to Eq. ( p3.12 ) and D'^" ^pip) is the most general 
solution to the corresponding homogeneous equation. Correspondingly [see Eq. ( 5.1l| )], we can write G[^""'(p) — 



Gp^^^ip) + GJj {p). To find the particular solution, we try an ansatz of the form 

D;\p{p) = diip) <5f„5^) + d^{p)p'P^''r,^p. 
Substituting this ansatz into Eqs. ( ^.12| ), it is easy to see that it solves these equations if we take 



di{p) = 



Flip) 



d2{p) = - 



F2{P) 



Fi{p)F^{p) 



with 



F3(p) = F^{p) + 3p^F2{p) = 1 - 48^Gp2 H^{p) _ 2^ 



(5.13) 



(5.14) 



(5.15) 



and where the notation [ ]r means that the zeros of the denominators are regulated with appropriate prescriptions 
in such a way that di{p) and d2{p) are well defined Lorentz invariant scalar distributions. This yields a particular 
solution to the Einstein-Langevin equations: 



G^;^^-'{p) = 16nGD;''^^{p)C^{p), 



(5.16) 



which, since the stochastic source is conserved, satisfies the conservation condition. Note that, in the case of a 
massless scalar field, m = 0, the above solution has a functional form analogous to that of the solutions of linearized 
semiclassical gravity found in the Appendix of Ref. [^. Notice also that, for a massless conformally coupled field, 
m = and A^ = 0, the second term in the rig ht ha nd side of Eq. (5.13) will not contribute in the correlation functions 



(5.10), since, as we have pointed out in Sec. IV B, in this case the stochastic source is "traceless." 



Next, we can work out the general form for ^isip)' which is a linear combination of terms consisting of a Lorentz 



invariant scalar distribution times one of the products (5^^(5^^, p^P^^rj, 



7a/3, Tl'^^Vo^p, V^-'p'Po.p, p'P;' and p-'P'^'^p-'P, 



a/3- 
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However, taking into account that the stochastic source is conserved, we can omit some terms in D'^ ^^{p) and simply 
write 



with 



(5.17) 



(5.18) 



where hi(j>), h2{p) and h^i^p) are Lorentz invariant scalar distributions. From the fact that D'^ ^^{p) must satisfy the 
homogeneous equation corresponding to Eq. (5.12), we find that hi{p) and /i3(p) have support on the set of points 
{p^} for which Fi{p) = 0, and that /i2(p) has support on the set of points {p^} for which Fi{p) = or F3(p) = 0. 
Moreover, the conservation condition for G''^^^^{p) implies that the term with h^{p) is only allowed in the case of a 



(5.19) 



massless conformally coupled field, m — O and A^ = 0. From ( |5.7| ), we get 

(G$:'^"^(p)r^(p')>c = {2nflQnG5\p + p')D>;:'(p)N'"^'^P{p). 



Note, from expressions (5^) and (5^), that the support of N^^'^"^{p) is on the set of points {p^} for w hic h —p ^ > 
when m — 0, and for which —p^ — Ani? > when m 7^ 0. At such points, using expressions dStel), ( |5l5| ), ( pll ) and 
(4.3), it is easy to see that Fi{p) is always different from zero, and that F3(p) is also always different from zero, except 
for some particular values of and /3: 

a) when m — 0, A^ = and /3>0; 

b) when m^O, 0<A^<(1/12) and P^{A^/32Tr^)[n/{Gm^) + 1/36]. 

In the case a), F^lp) — for the set of points {p'^} satisfying — — l/(967rG/3); in the case b), F^^p) = for 
{p'^} such that —p^ — m^/(3A^). Hence, except for the above cases a) and b), the intersection of the sup- 
ports of N'^'^°'^{p) and D^^ y^^{p) is an empty set and, thus, the correlation function ( 5.19| ) is zero. In the 
cases a) and b), we can have a contribution to (5.19) coming from the term with /i2(p) in ( 5.18|) of the form 
D^" p^{p)NP''°'f^{p) = H3{p;{C})p'^P'"' N°''^Pp{p), where H3{p;{G}) is the most general Lorentz invariant distribu- 
tion satisfying F-i{p) H-i{p; {C}) = 0, which depends on a set of arbitrary parameters represented as {C}. However, 
from (^), we see that N°'^Pp{p) is proportional to 0(— p^— 4m^) (l + 4m^/p^)(^/^) (3A^ + to^/p^)^. Thus, in the case a), 
we have N°'^''p{p) — and, in the case b), the intersection of the supports of N^^^pip) and of H^{p] {C}) is an empty 
set. Therefore, from the above analysis, we conclude that &^^^'^ {p) gives no contribution to the correlation functions 

( pO| ), since (Gi^''^''(p) r'^(p')lc = 0, and we have simply ^^^''"^(a;, a;') = (G<,'"'''(x)G<,i'"^(a;'))c, where G^^^^^ix) is the 
inverse Fourier transform of (^.1(: ) . 

The correlation functions ( |5.10 ) can then be computed from 

{&^''"'{p)&r^{p))c = 64{2nfG'S\p + p')D^''p^{p)Df,^{-p)N''-^^p). (5.20) 

It is easy to see from the above analysis that the prescriptions [ ]r in the factors Dp are irrelevant in the last expression 
and, thus , they can be suppressed. Taking into account that Fi{—p) = F*{p), with Z = l,2,3, we get from Eqs. ( 5.12 ) 
and dsl^ 

SHp + p') 



(G^^'P^ip) &^"^^\p'))c = 64 (2^)6 G^ 



\Fi{pW 



F3[P) 



p^p-^m-P (p) + \IlM}^p^pt^-p^p-fiNP, 



F^'ip) 



FzipW 



ip) 



(5.21) 



This last expression is well defined as a bi-distribution and can be easily evaluated using Eq. (5.8). We find 

2 



(G<^)''''(p)G<^>"''(p')), 



45 



(2.)^G^ii^e(V- 



\Fiip)\ 



-Am 




a0\ 



10 3 AC +^ (p2)2pM.pa/3 



1 - 3p' 



F2{P) 



Fsip) 



(5.22) 
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To derive the correlation functions (5.1C), we have to take the inverse Fourier transform of the above result. We 
finally obtain 



(5.23) 



with 



m 



gA{p) = 0i-p'-Am')Jl + A— 1 + 4 — 



1 



\Fi{p)\ 

2X2 1 



l-3p 



2 F2{p) 



(5.24) 



and JP^''"'^ = ajP^^^JPi?^" - T^^^'T^'^, and where i^;(p), Z = 1,2,3, are given in (IJ) and (|1|). Notice that, for a 
massless field (m = 0), we have 



Flip) ^l + l6nGp^Hip;fl^), 
16 r 

F2ip) = -—7TG[il + 180 AC') i7(p; /I^) „ 6T 
Fsip) = 1 - 487rGp2 [gO AC^iJ(p; /i^) - 2T 



(5.25) 



with fi = fj, exp(19207r^a) and T = /3 — GOA^^a, and where H{p]iJ?) is the Fourier transform of H{x;iJ?) given in 



B. Conformal field case 



The above correlation functions become simpler when the scalar field is massless and conformally coupled, i.e. 



when m — Q and A^ = 0, since in this case C/b(x) = and Qa(j>) reduces to Qa{p) 



-p )\Fi{p)\ . Introducing the 



function ip{x] A) = [1 — xln(Ax/e)] + X ^ with x > and A > 0, 5a (2;) can be written as 



Ga{.x) 



(1207r) 



3/2 



1 r°° 



q| sm 



Vl207r 
Lp 



|x||q| 



dg" cos 



Vl207r 



^(-g2;A)' 



(5.26) 



where Lp = vG is the Planck length, A = 1207re/(Lp/i2), and we use the notation x^^ — (a;°,x) and = (<z°,q). 
Notice that, if we assume that jl < Lp^, then A ^ 10^. For those values of the parameter A (and also for smaller 
values), the function ip{x'. A) has a minimum at some value of x that we denote as xo(A). This can be found by solving 
the equation tt^xo = [1 ~ Xoln(Axo/e)] [1 + ln(Axo/e)] numerically (discarding a solution XAf(A) < Xo(A), at which 
the function ip{x', A) has a maximum). Since the main contribution to the integral ( 5.26 ) come from the values of —q^ 



around —q^ = xo(A), fixi A) can be approximately replaced in this integral by (papixi ^) = ~ '^(•^) x]^ + ^'^X^ = 
[k2(A) +7r2]x^ - 2k(A)x+ 1, with k{X) = ln(Axo(A)/e). For (A/5) - 10^ - 10^ we have k ~ 10. 
Let the spacetime points x and x' be different and spacelike separated. In this case, we can choose an inertial 



coordinate system for which [x 
written as 



(0,x — x') and G'^'^°'^{x,x') will be a function of x — x' only that can be 



(x - x') - gr'^^i^ - x') + gr^'^ix - 



(5.27) 



with 



c;r"''(x)^^G2^^/"''/,(x), 



1,2,3, where /i(x) = gAix)],,^^^^..^ ^i^) = {dlfQA{x) 



la;f' = (0,x) 



h{-^)^{dlYgA{x)\ 



a;A' = (0,x) ■ 



(5.28) 
and jF^=^"/5 are 



some differential operators. Note that the terms containing an odd number of derivatives are zero. The differential 
operators Tl"^""^ are given by jr^^''"^ = sp^Cap/S).. _ p^i/pa/s^ ^j^j^ -j^t.v ^ (^M^jy _ Sf'^S'^i) d^dj. The non-nuU 

components of the remaining operators are = 39*9^ - (5*^A, = \ [d'd^ + 35*JA), T^^''^ = -5'^ 5^^ + 3(5'('=(5')J, 
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jq^''^ = 2 ((5'^^'=' - 3(5*('=(S')J) A -S^^d^d^ - S'^^d'd^ + 3 {S^^^d^^d^ + where AEE^'^aiSj is the usual (EucUdean 

space) Laplace operator. From the above expressions, we can see that t/'^°*'*(x — x') — Q'^^^^{yi — x') = 0, but the 
remaining correlation functions Q>^'^°'^{x — x') are in principle non-null. 

With the approximation described above, the integrals /a(x) can be written as 



(-1)°+! /120^Y 1 



d|q| sin 



|x||q| 



|q|^a(|q|), 



(5.29) 



where 



^i(|q|) 

^3(|q|) 



Vap(-g^;A)^ 



dq" 



K2(A)+7r2 



|q| 



lql '^ap(-g2;A)' 
1 



(5.30) 



Noting that ips_p{—q'^; A) has four zeros in the complex g° plane at ±p(|q|), ±p*(|q|), where p{s) (we make s= |q|) is 
the complex function with 



Rep(s) 
Imp(s) 



[(k2 + 7r2) s2 + + 7r2 ± (k2 + ^2) _^ ^ 
2(^2 +7r2) 



1/2 



we can decompose 
1 



V3ap(-g2;A) 4(^2 + 7r2) |p|2 Rep 



g° + 2Rep 



((j° - 2Rep) 



(g0)2 + 2RepgO + \p\^ (g0)2 - 2RepqO + 



(5.31) 



(5.32) 



and then we can perform the integrals Ja{s), a = 1, 2, 3. The results for these integrals can be found in Appendix 

Next, to carr y on w ith the calculation, we need to introduce some suitable approximations for the functions Ja(s) 
in the integrals ( 5.29 ). In order to do so, we study the behavior of these functions for small and large values of s. For 
s Ji{s), we find that it can be well approximated by an arctan function. In fact, on the one hand, s Ji(s) tends very 
quickly to a constant limiting value lims^oo s Ji(s) — a/4, where a = 1 + (2/7r) arctan(K/7r). On the other hand, for 

small values of s, we havesJi(s)~ [Vl207ra/(27r6)] s + 0{s^), with = {'ia/n^)[l5TT{V n'^ + n'^ - k)] . Hence, we 
can approximate 



s Ji(s) 



2tt 



arctan 



VUOn 



(5.33) 



Performing the integral /i(x) [see Eq. (5.29)] with this approximation, we get, for |x|7^0, 

15 a 1 
T 



^i(x) 



-b\x\/Lp 



(5.34) 



The function J2(s) behaves as J2(s) — (a/4) s + 0{s ^ Ins) for large values of s, and as J2(s) ~ (a/4) (1207r) ^/2 ^ _|_ 
0(s2), with 7 = 240 (k^ + 7r2)~i/2 ^-i^ fQj. gmall values of s. This function can be well approximated by 



Ms) - i 



1207r 



1/2 



and, substituting the last expression in the integral /2(x) [see ( 5.29| )], we obtain, for [xIt^O, 

15 a 72 



/2(x) 



^2 Lj, |x|2 



if2(7|x|/Lp), 



(5.35) 



(5.36) 



where K^{z) denote the modified Bessel functions of the second kind. For J3(s), we find that Jsis) ~ 
(a/4)s3 + ©(sins) for large values of s, and that Jsis) ~ (a/4) (1207r)-3/2 + 0(s), with d = A{k^ + 

7^2)^1/2 |^4507r6^^ (2k — \/K^~fn^)] , for s small. With the approximation 
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3/2 



Ms) 



1207r 



we can compute the integral /3(x) [see ( 5.29| )] for |x| 7^0, and we find 



(5.37) 



(5.38) 



Numerical calculations confirm that the above approximations are reasonable. For k ~ f 0, we have a, 6, (5 ~ 1 and 
7 - fO. 

The resuhs ( jsl^ ), (|53|) and (|38|) are now rea dy to be substituted into ( ^.28[ ), from where we can compute the 
different contributions to the correlation functions ( 5.27|) . Using the relation [d/ dz)'yz^" Ki,{z)'\ =—z~'^K,y+i{z), and 
defining df, = 6 |x|/_Lp, cr^=7|x|/Lp, as = S Lp, we get, after a rather long but straightforward calculation, the 
following results for the non-zero components of tj^''"''(x) [with [xIt^O]: 



c?r°(x) 



2 e- 



3tt Lj, al 

1 ab^ e-"^ 



^r(x)=-^c^^(x), 



12 


24 


24" 




+ — 


+ — 


5 

h 


16 

— + 


32 









7 24 48 48 
—+—+—+— 

C^b f^b cr;j 



1 ab^ e-"^ 



CTb (T^- CT^3 (Tj, 



ct;^ (t;? ct? 



|X|^ V / V 

2 ^ , k if. 14 87 279 384' 
■ x'x^x^x'' (IH h — + ^ + ^ 



6 18 48 
(^l 4 4 



^2 V 7-4 



1 07^ i4r4(cr^) /^^ x^x^ 



(Jb (Jb 



1 076 i4r4(cr^) /^a;'a;J 
1 076 i^4(cr.^) 



^ 07^ 



37r Lp 



2 (^S'^6''^-3S'^^S^^^^- (s'^x''x^+S''^x'x'-3S''^''x^')x^ 



-3S^^''x'\' 



4 ^^^3(a^|'^.^.._3^.(.^0.A 
37r Lp at V / 



(5.39) 



TT Lp '^S ' 

Note that, for a~:^\, we have the following asymptotic expansions for the modified Bessel functions in the above 
expressions: 



. 63 1 ( \ 



i^3(a)^,/^e-^ 



1 35 1 f \ 



(5.40) 



C. Correlation functions for the metric perturbations 



Starting from the solutions found for the linearized Einstei n ten sor, which are characterized by the two-point 
correlation functions ( 5.23| ) [or, in terms of Fourier transforms, ( 5.22| )], we can now solve the equations for the metric 
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perturbations. Working in the harmonic gauge, d„h^^^ = Q (this zero must be understood in the same statistical sense as 
above), where h^^ = hf^,j—{\/2)r]^^h, and using Eqs. ( 2.19 ) and (El), these equations reduce to □/i'"'(a;) = — 2G''^'^''(a:), 



or, in terms of Fourier transforms, p^hi^'^{p) = 2(5*^'^" (p). As above, we can write /i'^" = {h'^" 
solution to these equations with zero average, and the two-point correlation functions are given by 



n 



ix,x') = {h^^-{x)h^f\x')),-{h^^-{x)),{h"f'{x')), = m^ixyhfix')) 



where hf^ is a 



(5.41) 



We can now seek solutions of the form hf'^ (p) — 2D{p)G\''^ '^'^ (p) , where D{p) is a Lorentz invariant scalar distribution 
in Minkowski spacetime, which is the most general solution of p^D{p) = 1. Note that, since the linearized Einstein 
tensor is conserved, solutions of this form automatically satisfy the harmonic gauge condition. As above, we can 
write D(ji) — [l/p'^]r+Dh{p), where Dh{p) is the most general solution to the associated homogeneous equation and, 
correspondingly, we have h'^'^ (p) — hp'^ (p) + h'^^{p). However, since Dh{p) has support on the set of points for which 
p2 = 0, it is easy to sec from Eq. ( ^.22[ ) [from the factor e{-p^ - -im^)] that (%''(p)G^'>"'^(p'))c = and, thus, the 

two-point correlation functions (5.41) can be computed from {h''^'' {p)hf {p'))^^ {h'^'' {p)hf {p'))c. From Eq. ( 15^2^ ) 
and due to the factor 9{—p^ — 4m'^), it is also easy to see that the prescription [ ]r is irrelevant in this correlation 
function and we obtain 



{-h';;-'{p)-hf{p')) 



{G^;'^''{p)G^;^'^f'{p')) 



(5.42) 



where {G'p''^''{p) G<,^'"^(p'))c is given in (5.22). The right hand side of this equation is a well defined bi-distribution, 
at least for m 7^ (the 9 function provides the suitable cutoff). In the massless field case, since the noise kernel is 
obtained as the limit m — > of the noise kernel for a massive field, it seems that the natural prescription to avoid 
the divergencies on the lightcone = is a Hadamard finite part (see Refs. for its definition). Taking this 
prescription, we also get a well defined bi-distribution for the massless limit of the last expression. Finally, we find 
the result 



4.77- 

n^''°'f^{x,x') = — G'^ T^^"'^^ nA(x 
45 



^^G^T^'^Tfn^ix-x'l 



(5.43) 



where ?^a(p) = [iKp^f] Qa{p) and Hb(p) = [l/{p^?] Gsip), with GAip) and geip) given by The two-point 

correlation functions for the metric perturbations can be easily obtained using hfj_^ = h^^ — (l/2)'i]^i,h". 



VI. DISCUSSION 



Our main results for the correlation functions are ( 5.23 ) and ( 5.45 ). In the case of a conformal field, th e cor relation 
functions of the linearized Einstein tensor have been explicitly evaluated and the results are given in (5.39). From 
the exponential factors e^*^ in these results, we see that the correlation functions of the linearized Einstein tensor 
are in this case characterized by correlation lengths of the order of the Planck length. A similar behavior is expected 
for the correlation functions of the metric perturbations. Hence, as expected in this case, the correlation functions 
are negligibly small for points separated by distances large compared to the Planck length. At such scales, the 
dynamics of gravitational perturbations of Minkowski spacetime can be simply described by semiclassical gravity 
31 4^,^,^,^. Deviations from semiclassical gravity are only important for points separated by Planckian or 



sub-Planckian scales. However, for such scales, our results (5.39) are not reliable, since we expect that gravitational 
fluctuations of genuine quantum nature to be relevant and, thus, the classical description breaks down. It is interesting 
to note, however, that these results for correlation functions are non-analytic in their characteristic correlation lengths. 
This kind of non-analytic behavior is actually quite typical of the solutions of Langevin-type equations with dissipative 
terms. An example in the context of a reduced version of the semiclassical Einstein-Langevin equation is given in 
Ref. [|0|). 

For background solutions of semiclassical gravity with other scales present apart from the Planck scales (for instance, 
for matter fields in a thermal state), stress-energy fluctuations may be important at larger scales. For such back- 
grounds, stochastic semiclassical gravity might predict correlation functions with characteristic correlation lengths 
much larger than the Planck scales, so as to be relevant and reliable on a certain range of scales. It seems quite 
plausible, nevertheless, that these correlation functions would remain non-analytic in their characteristic correlation 
lengths. This would imply that these correlation functions could not be obtained from a calculation involving a 
perturbative expansion in the characteristic correlation lengths. In particular, if these correlation lengths are propor- 
tional to the Planck constant h, the gravitational correlation functions could not be obtained from an expansion in h. 
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Hence, stochastic semiclassical gravity might predict a behavior for gravitational correlation functions different from 
that of the analogous functions in perturbative quantum gravity . This is not necessarily inconsistent with having 
neglected action terms of higher order in Ti when considering semiclassical gravity as an effective theory ||2^ . 

We conclude this section with some comments about a technical point on the obtained solutions of stochastic 
semiclassical gravity. It concerns the issue that the Einstein-Langevin equations, as well as the semiclassical Einstein 
equations, contain derivatives of order higher than two. Due to this fact, these equations can have some "pathological" 
solutions {e.g., "runaway" solutions) which are presumably unphysical | |45| , ^ , ^ , p5| . Thus, one needs to apply some 
criterion to discern the "physical" from the unphysical solutions. However, as it is discussed in Ref. ||25| (see also 
Refs. 1^^), even in the context of "pure" (non-stochastic) semiclassical gravity, this is still an open problem. Two 
main proposals, both based in the works by Simon have been made concerning this issue: the "perturbative 

expandability" (in Ti) criterion |4^,^,^ and the "reduction of order" procedure . 

The first proposal consists in identifying a subclass of "physical" solutions which are analytic in the Planck constant 
h. This proposal has been successful in eliminating the instability of Minkowski spacetime found by Horowit z pO| , pT| . 
However, on the one hand, this proposal seems to be too restrictive since, as it has been pointed out in Ref. \2b\, one 
could not describe effects such as the continuous mass loss of a black hole due to Hawking radiation. On the other 
hand, there can be situations in which the formal series obtained when seeking approximate perturbative solutions 
(to a finite order in K) does not converge to a solution to the semiclassical equations j2^. In our case, if we had tried 



to find solutions to Eq. (5^) as a Taylor expansion in %, we would have obtained a series for G^^l{p) which, as the 

above solutions, would be linear and local i n ^a0 {p), but whose corresponding two-point correlation functions for the 
conformal field case would not converge to ( 5.23| ). 

The "reduction of order" procedure provides in some cases a reasonable wa y to modify the semiclassical equations in 
order to eliminate spurious solutions. But, as it has been emphasized in Ref. (25|, it is not clear at all whether a reduc- 
tion of order procedure can always be applied to the semiclassical Einstein equation (and how this procedure should be 
applied). For the Einstein-Langevin equation, this issue has not been, to our knowledge, properly addressed. A naive 
application of the prescription to Eq. ( |5.2| ) seems to downplay the role of the dissipative terms with respect to the noise 
source. In fact, to lowest order, we obtain G*^'^'' = IGttG^^'^, where there is no contribution of the dissipation kernel. 
From this equation, we get the well-known result (G"''"'')c = ||2|,||, and also g^'''°'f\x,x')^{lQ^:f L%Ni"'"l^ {x,x'). 

and (|t|), this gives g^"'°'l^ {x,x') = {2llb){L%/-K^) \{l/Q)T^''°'^ + 



For a massless field, using Eqs. 



Vf 



X 



1/{{X 

n'"'"'^ix,x 



( P-ioD 

v2 



(4.9) 



^3x4 



6%X 



the harmonic gauge 



') ^ {ATT/A5)L%Tt^''"'^JA{x - J 



For the two-point correlation functions (5.41), we get, in 

+ (327r/9)L| 
4m2) 



T^^" Tf^ Iji{x - x'), with lx{p) = 



p" ^im" ) [p" ) " L+im-^/p'" [^L +4:771" /p"j and J-B(p)^t^(— p" — 4m") (v")' 
paring the last results for the massless case with the ones obtained in Sec. [V|, we note that the main qualitative feature 
is the absence of the exponential factors e~'^ , which make the two-point correlation functions to decay much more 
slowly with the distance, i.e., like a power instead of an exponential law. This fact is due to the lack of dissipative 
terms in the reduced order equations. The conclusion is that one should probably implement a more sophisticated 
version of the reduction of order procedure so as to keep some contribution of the dissipation kernel in the reduced 
order equations. 

For these reasons, in our wo rk w e have not attempted any of these procedures and we have simply sought some 
solutions to the full equations (5^). Our solutions for the conformal field case have the physically reasonable feature 
of having negligible two-point functions for points separated by scales larger than the Planck length. 
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APPENDIX A: THE KERNELS FOR A VACUUM STATE 



The kernels for a vacuum state can be computed in terms of the Wightman and Feynman functions defined in 
Eq. (2.6) using 



I it{x) i^f{y) |0) = 4 {Nf^''{x, y) + iHf^^^x, y)) = V^G+(x, y) v'G+(x, y) 
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where V is the differential operator 



and 



+ 2I?'"'(x^^G,„(x,2/) ^''G,^{x,y)) + 2Vf{yJ^G,,^{x,y) X^^G^Jrr, y)) 
+2I?"''I?;^(G,2(x,y)) + i [g'"'(x)(^^G.„(a;,y)v/+\^'^G.„(:r,y). 



(Al) 



(A2) 



(A3) 



APPENDIX B: MOMENTUM INTEGRALS 



Some useful expressions for the momentum integrals in dimensional regularization defined in (3.14) and ( 3.15| ) are 



(47r)2 
0, 



47r^2 



^'1-2 



(47r)2 (n-2) 



K„ + 0(n-4), 



TP- 

^0,1 ~ m ri ^ 



Jn{p) = [2k„ + 0(p2) + 0(n-4) 



2 

4 

Jn{p) 



1 + 4- 







/ (n-l)J 


2 p2 



p^'p 



n-1 



4 p2 



3p>V 



pi^p'^p^pP - ( 1 + 4 



(n-1) 



02 

, , , ^—r (P'"'p" + P'^^'p" + P°"'p'') 
p2 J [n—lj ^ ' 



p2 y (n — 1) 



{P^'^p^p^ + P'^^'p^'p^ + P'^f^p^'p"' 



P^yp/^ + Pi'^'p-'p"' + P^f^pfp") + ( 1 + 4 



2\2 /^2\2 



{P'f 



(n2_l) 



V pQ/3 



8 p2 



pt^p^pO'pli + 



71 



1 4 m" 



n — 1 n p 



p2 ^ p-^^pMp/? 



(Bl) 
(B2) 
(B3) 
(B4) 

(B5) 
(B6) 

(B7) 
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p^^pfp" + P''"p''p>^ + pff^p^p" + p^f^pfp") 



1 4(2n-l) 



n^ — 1 71(71^ — 1) p^ 



where p^P'"' = t?'"'^^ _ ^t^^u ^ defined in ( pl6| ), 

(j){p^) = da ln^l + ^a(l-a)-iej = -ivr ^(-p^ - 4m2) ^1 + 4 ^ + (y3(/), 
with e-^0+, and 



da In 



1 + ATa(l-a) 



'1 + 4 In 



i + 42i; + i 



1 + 4 ™l _ 1 



1 + 4- 



-1-4- 



We can also write (jjijP') in a more compact way as 



</.(p2) = ^2 + Jl + 4— hi 



^1 + A{m? - ie)lp^ + 1 
^/l + 4 (m^ - it) Ip^ - 1 



Other useful integrals in momentum space defined in (3.5) are 
1 ' 



4(27r)3 



(V)^(-/-4™^)Wl + 4^, 



m 
4 



( 1+4— 



8 

m 

16 



2 \ 2 

- ( 1 + 4 ^ j ^ (F^'^p" + P^"p'' + P" V 



2 \ ^2 



pf^p'^pO'pP „ I 1 + 4 ^ j ^ (Pt'^p'^pf^ + P'^'p'^pf^ + P''^p''p°' + P'^'p'^p'^ 



pt^0p<^pa pa 13 11 



P'-p- 



2 \ 2 / 2\2 

1-1-4 I ' ( pt^" P°'P -f pA*" P''/3 -I- pM/3 pi^ 

p"^ J 15 



(B8) 



(B9) 



(BIO) 
(Bll) 

(B12) 

(B13) 
(B14) 
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APPENDIX C: PRODUCTS OF WIGHTMAN FUNCTIONS 



For the products of derivatives of Wightman functions involved in the calculations of Sec. Ill A, we obtain the 
following expressions: 



d'^p 

e 



4, 



A-(.x)--(2.)^^(2^), 
5''A+(x)a^A+(x) = (27r)2y'^^^^^ 

9''a^A+(a;)5"<9''A+(:E) = -(2^)^ / 

J v-T^r 



dp e-'f-[/^(p)p'^-/^''(p)], 



-ipx 



I'"'ip)p°'p^ - 2/'"'("(p)/) + If'^^^ip) 



(CI) 
(C2) 

(C3) 
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with I{p), I^'(p), I'"'{p), I^"'°'{p) and I^"'°''^{p) given by Eqs. ( p312| )-( p316[ ). From these expressions, using the results 
of Appendix H, we obtain 



^ !(„) ^'nljTt^'^ f jUl 

(27r)4 3 " 7 (27r)4 



1 + 4^ /(p), 



Hp) 



12 

^2 



d p / ^ . 



(27r)4 



Hp) 



~2 r ^4 

60 ^ ^ ^ V (27r)4 



1 + 4- 



1 + 4-^) /(p). 



(C4) 
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APPENDIX D: PRODUCTS OF FEYNMAN FUNCTIONS 



For the products of derivatives of Feynman functions that we need for the calculations of Sec. [II B , we obtain the 
following results: 



(27r)" 

-4)5^A,„(x)a^A,„(x) = - 



(2^)^ 



(27r)« 



-4)A,„(0) = -/o„, 



(27r)" 



e^^-(//o„+2pX.+CA<) 



Using the results of Appendix O, we find from the above expressions 
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15 \ n J 



i 1 



^ {p^ + 6m^) {P>"'p°'p^ + P°'l^p>'p'') 



(47r)2 9 



i 1 



(2 (p2)2 + 20mV + 45m'') (pt'^P"'^ + 2pt'(°'pl^>) I + 0(n - 4), 
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where P^" is the projector orthogonal to p^ defined above 



APPENDIX E: LINEARIZED TENSORS AROUND FLAT SPACETIME 

Some curvature tensors hnearized around flat spacetime are given by the foUowing expressions: 



with 



Q/3 



and 



In four spacetime dimensions, the hnearized Weyl tensor is given by 

- rj^'I^T^P^Pd^d" - Tj''f^r]°"'dPdP - f]"'' rjP"^ d^" - f]'"'r]PpTj^''a + if'^r^PPd'^d" 
+ 7j''°'7j'^''dPdP - T^P^ifd'^d'^ - r^Pf^r^^PifU + r^t^f^ t^'' p d"^ d" + j^pI^ ri"" d" dP) 



APPENDIX F: THE INTEGRALS Ja{s) 



For the integrals Ja(s), a=l,2.3, defined in (5.30), we find the following results 
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. I r 2 11/2 

where p is a function of s given by expressions (5.31), which give 
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